Spin accumulation with spin-orbit interaction 
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Spin accumulation is a crucial but imprecise concept in spintronics. In metal-based spintronics 
it is characterized in terms of semiclassical distribution functions. In semiconductors with a strong 
spin-orbit coupling the spin accumulation is interpreted as a superposition of coherent eigenstates. 
Both views can be reconciled by taking into account the electron-electron interaction: a sufficiently 
strong self-consistent exchange field reduces a spin accumulation to a chemical potential difference 
between the two spin bands even in the presence of spin-orbit coupling. We demonstrate the idea 
on a clean two-dimensional electron gas (2DEG) by showing how the exchange held protects a spin 
accumulation from dephasing and introduces an easy-plane anisotropy. 
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Metal-based spintronics [T] has evolved into a mature 
field in which spin phenomena are routinely exploited in 
versatile applications [2]. However, integration of spin- 
based functionalities into semiconductor circuits is still 
a pressing challenge. Much of recent research in this 
area has been motivated by device concepts, such as the 
seminal Datta-Das transistor [3], which requires injec- 
tion and detection of spins by ferromagnetic contacts to 
a narrow channel of a two-dimensional electron gas with 
gate-controlled spin-orbit interaction (SOI). In spite of 
progress to inject, modulate, transport, and detect spin 
polarization ail-electrically [4 13, 0, S B 0] well as 
evidence that the SOI can indeed be tuned by exter- 
nal gates p^, the route to a working spin transistor 
appears to be still full of obstacles. In the meantime, 
many insights have been obtained on the spin accumu- 
lation and its dynamics by optical methods, especially 
time-dependent Kerr and Faraday rotation spectroscopy 
[nl, [H, Q [H, Large Rashba splitting have 

been observed at metal surfaces by angle-resolved pho- 
toemission [17] , which attracts a lot of attention recently 



We define a spin accumulation as a non-equilibrium 
spin-polarized state injected optically or electrically into 
a non-magnetic material. In metal-based spintronics a 
spin accumulation is synonymous to a chemical poten- 
tial difference between spin up and down bands [20]. 
However, in semiconductors the SOI prominently affects 
the electronic structure and transport properties |2l|. 
A spin accumulation is then interpreted as an intrinsi- 
cally time-dependent quantum superposition of coherent 
eigenstates. This difference is not just a semantic ques- 
tion but essential for the functionality of spintronic de- 
vices. Here we offer a unified mean-field theory for the 
spin accumulation in both metals and semiconductors. 
Spin can be injected either slowly, e.g. by a ferromag- 
netic contact with small electric bias, or rapidly, e.g. by 
pulsed optically induced excitation. We start below with 
a description of spin-accumulation eigenstates that are 



accessible by adiabatic excitation followed by a brief dis- 
ucssion of the spin accumulation dynamics of rapidly ex- 
cited states. We illustrate the general ideas at the hand 
of a 2DEGs with Rashba SOI [22|, in which the disorder- 
scattering lifetime broadening is much smaller than the 
spin-orbit splitting at the Fermi-level. 

Let us consider an infinitely extended homogeneous 
2DEG. To leading order in the electron wave vector 



k = —iV the Hamiltonian including the SOI is |21| 
Hq — fi^/c^/2m* + a{axky — Cykx)^ where m* = mrrrie 
is the effective electron mass, mg is the bare electron 
mass, ai are the Pauli matrices, and a is the Rashba SOI 
strength parameter 0]. Electron-electron interactions 
are treated within the density-matrix functional theory 
(DMFT) [23] which is a generalization of the Hohenberg- 
Kohn-Sham density-functional theory 0] that can han- 
dle excited states. Compared to the Hartree-Fock (HE) 
method, exchange and correlation effects in the DFMT 
can be treated within local approximations. The re- 
duced density matrix is r(z,z') = Xli^ ^*Xi (z) xl (^O^ 
where Xi eigenstates of the Kohn-Sham Hamilto- 
nian (natural orbitals), < < 1 are the corre- 
sponding eigenvalues (natural occupation numbers), and 
z = (r, cr) is space-spin coordinate. We define Fg as 
the subset of all density matrices which correspond to 
a given electron density p (r) = ^ixl (z) Xi (z) and 
spin-polarization s (r) = ^ixl (z) crXi (z) /N, where 
N is the total number of electrons. The density-matrix 
functional is defined via minimization of the total en- 
ergy E[Ts] = min^[r^](^[Fs]|i^|^[Fs]) in the space of ah 
many-body wave functions that correspond to a given Fg. 

We now assume that the exact density matrix can be 
generated by a non-interacting system of pseudo particles 



[Ho + Vea:t -^Vh^ Kc]^i = 



(1) 



where Vext^ Vr^ and V^c are the external, the Hartree, 
and the exchange-correlation potential, respectively, such 

that Fs(z,zO = Ef /^^^(Z)4(Z0 with /, = {0,1} 

and l/xc(z,z', [Fs]) = ^£^xc[rs]/^Fs(z, z'), where E^^ is 
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the exchange-correlation energy. HF-calculations for 
the Rashba Hamiltonian, fohowing Ref. [1^, confirm 
that the effect of the SOI on the exchange potential 
is negligible for small spin polarizations. With the lo- 
cal approximation, we finally arrive at l^c(z,z', [Fg]) ~ 
^(r — r'){Vo{p^ s) + Jxc(p, s)s • a), where s = ss, and 
J^c{p-,s) J{p)s for small s) is the modulus of the 
exchange-correlation field vector. The scalar Vb can be 
dropped in homogeneous systems. We may approximate 
J by the HF exchange energy of the strictly 2DEG [25| 
Kcc ^ J^{p-) ' ^ J{p)^ ' ^ = — V^mr/ (7rrse^)Ry s • a, 
where Vg = mr/{aB^^/^) is the dimensionless density 
parameter, e is the relative dielectric constant of the 
medium, Ry = 13.6 eV, and as = 0.53 A. The effec- 
tive Hamiltonian is then 



H (s) = Ho^ Js- a, 



(2) 



where J < is the effective exchange potential, in which 
correlations can be included using published parameter- 
izations of the correlation energy for a non-SO coupled 
2DEG [27|] . For typical electron densities p the exchange 
energies (a few meV) are of the same order of magni- 
tude as SOI energies at the Fermi-level in III-V and 
II-VI ^28] semiconductor-based 2DEG's. 

The eigenstates of the Hamiltonian in the non- 
polarized ground state are split into two bands with a 
chiral spin pattern (Fig. [T](a-b)). The exchange field de- 




FIG. 1: (Left panel) Spins (a) and energies (b) of the ground- 
state spin-split bands in a non- interacting Rashba 2DEG. The 
inner and outer circles correspond to the Fermi energies of the 
spin-split bands. (Right panel) Modulation of the electronic 
structure in the presence of an in-plane spin accumulation 
s by its the exchange field. Circles in (c) are fixed-energy 
countours. The shifted occupation number distributions that 
minimize the energy are shown in (d). 

forms the electronic bands and spinors as shown in Fig. 



[TJc) for in-plane and in Fig. [2ja-b) for perpendicular di- 
rection of an injected spin accumulation. SO split bands 
of the surface states of in-plane magnetized Gd films have 
been found to be deformed by the exchange potential fi^ 
similar to Fig. [D^c). Our task is to find the self-consistent 
single-determinant eigenfunction of (s), which accord- 
ing to the DMFT is unique. Introducing the occupation 



spin polarized 
(0< sz « 1).^' 
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FIG. 2: (a-b) Spin direction of the lower spin-split (— ) band 
in the presence of a spin polarization normal to the 2DEG 
plane. The total spin polarization is determined by the shaded 
area between kp,- and - c) The ground state (solid line) 
and excited eigenstates (dashed line) with non-zero, small 
spin polarization perpendicular to the 2DEG surface are sep- 
arated with an energy gap. d) Band polarization of spin- 
accumulation states at fixed J — —2 meV. Material parame- 
ters are e 12.7, m* 0.067me, /? = 2 x 10^^/m^ 

numbers /kA = {0,1} of the spin-split states ^kA, the 
spin polarization reads 



S = ^/kA(0kAk|</.kA)/iV, 



(3) 



kA 



where A = {+, — } is the band index. The state we 
are looking for minimizes the energy under the con- 
straint ([3]), with the /kA as constrained variational vari- 
ables. We solve the problem either analytically in lim- 
iting cases or numerically by a stochastic minimization 
method, which uses penalty functions to fix the spin po- 
larization and the Metropolis sampling method to find 
the global energy minimum. An unequal occupation of 
spin bands, 7V_ = Xlk- /k- > ^+ can be parameterized 
by a chemical potential difference or a band polarization 
= {N- — N-^)/N > 0. Occupations can also shift in 
momentum space (Fig. Hl^d)). 

A spin accumulation in the plane of the 2DEG can 
be generated at minimized energy by shifting the Fermi 
circles, which induces currents [29] via a "spin-galvanic 
Hall effect": the minimum energy state at fixed Sx is 
associated with a charge current in the perpendicular y- 
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direction 



a / Jm* \ 



(4) 



Since electrons move with constant drift velocity, there 
are no intrinsic spin-Hall currents [30]. Non-equilibrium 
spin currents are induced, but vanish to first order in Sx- 
A spin accumulation perpendicular to the 2DEG sur- 
face can be generated by the exchange field that pops 
the in-plane spin textures out of the plane (Fig. El^b)). 
However, the SOI counteracts the spin- alignment and the 
band polarization must be increased from that of the 
ground state to support a spin accumulation [si^ . Conse- 
quently, excited eigenstates corresponding to a finite Sz 
are, in contrast to the in-plane case, separated from the 
non-polarized ground state by a finite energy gap 



Ea 



J 



(5) 



The divergence in at J ^ reflects the absence of an 
Sz component in the non-interacting Rashba model. The 
gap is shown in Fig. Wis^) in which Eq. ([5]) corresponds to 
the low a limit of the energy difference. This gap energy 
must be overcome to achieve spin-polarized eigenstates 
at arbitrarily small 7^ 0. Except for this singular be- 
havior at = we find that the energy of eigenstates 
is isotropic (to second order in s) to the numerical accu- 
racy for 5^ 7^ 0. We suspect that there may be a more 
general physical reason behind this out-of-plane isotropy. 
The contribution of spin-galvanic currents to E{s) is not 
significant for material parameters shown in Fig. O 

The maximum spin accumulation that can be accom- 
modated perpendicular to the 2DEG surface is deter- 
mined by the total spin polarization of a single occupied 
band (p5 = 1). In the exchange-only approximation the 
self-consistency criterion ([3]) can be fulfilled only when 



a<2\J\/kF = 2\J\/y 



0.32 eV nm/e, (6) 



Figure [21(d) shows the stability limit as a function of Sz- 
Per definition, eigenstates do not dephase. The dy- 
namics of the semiclassical spin accumulations discussed 
above is therefore governed by the Bloch equation s = 
-7s X Beff— s/Ti 0, in which Beff = —dE{s) /ds and 
Ti is the spin relaxation time. Due to the singular 
anisotropy of E (s) the Bloch equation is mathematically 
not well defined in the mean- field theory employed here. 
Qualitatively, the absence of an angle dependence of E{s) 
(for Sz 7^ 0) implies that exchange-stabilized spin accu- 
mulations do not feel an internal SOI field and precess 
only when an external magnetic field is applied. A spin 
accumulation exactly in the 2DEG plane is trapped and 
precesses around an in-plane external magnetic field that 
exceeds a threshold value governed by the energy gap ([5]) . 

A large phase space available for scattering processes 
makes a large spin accumulation susceptible to fast de- 
cay. Therefore the stability limit of eigenstates (j6|) is 



not a sharp phase boundary. For not too highly ex- 
cited systems, the Dyakonov-Perel [32] mechanism by 
random scattering at defects is believed to be the domi- 
nant source of finite Ti spin life times in clean systems. 
Since in an exchange stabilized 2DEG the precession in 
the SO field is suppressed, the efficiency of the Dyakonov- 
Perel mechanism is strongly diminished for systems in 
the clean limit considered here. The opposite (dirty) 
regime can be handled by spin-coherent kinetic H, 3^ 
and diffusion 



361 , 137j equations or numerical simula- 
tions [1^,^. 

The Datta-Das transistor is a spin valve consisting 
of a 2DEG spacer with transparent ferromagnetic con- 
tacts [3]. Even when the magnetizations of the two elec- 
trodes are parallel to each other, transport depends on 
the magnetization direction when exchange is taken into 
account. When magnetizations are oriented in the cur- 
rent direction, the spin accumulation can be injected 
into the 2DEG as an eigenstate with shifted distribu- 
tions (Fig. m^c)), which does not precess and, hence, 
does not react to a gate voltage that modulates the SOI. 
For magnetizations not in the 2DEG plane spin accumu- 
lation eigenstates are separated from the ground state 
by an energy gap and spin cannot be injected adiabati- 
cally at low energies. A non-adiabatic spin injection, on 
the other hand, could lead to a coherent superposition of 
eigenstates and spins would precess in the SOI field, as 
envisioned by Datta and Das. 

In optical pump and probe experiments, spin-polarized 
electrons and holes are generated by short resonant pulses 
of circularly polarized light, followed by fast thermaliza- 
tion and spin relaxation of the holes. A fast excitation 
creates a coherent superposition of individual spin eigen- 
states which dephases with time. We calculate the dy- 
namics of the spin accumulation from the initial state ipo 
for times t > from |V^(t)) = / e^^^^)^/^|V^o), where the 
Hamiltonian depends on t by the exchange field Js {t) . 
The state at t is solved iteratively 



JH{t-At)At/h 



(7) 

for short time steps At 1 fs. We assume that dephasing 
is fast compared to the spin relaxation processes so that 
the occupation numbers are unchanged. Figure [3] shows 
the time evolution of the spin accumulation, which has 
been excited at time t = into a coherent superposition 
of eigenstates. A single spin oscillates in the SO-field 
by a frequency ujso = 2ak/h. The spin polarization ex- 
cited over a finite band width is therefore expected to 
decay on the scale of the dephasing time T2, that de- 
creases with increasing s. However, a strong exchange 
field aligns spins along a common axis and synchronizes 
spin precession which protects the spin polarization from 
dephasing. The exchange-induced enhancement of T2 be- 
comes significant when the exchange-splitting, which is 
proportional to 5, becomes of the same order of magni- 
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FIG. 3: Oscillation and dephasing of a spin ensemble in a SO- 
field. The state is a coherent superposition of spin eigenstates, 
s{t = 0) = 10%, and a = 4 • 10"^^ eV m. 



tude as the spin-orbit splitting. Such an effect has been 
observed in experiments but can be explained by 
the exchange effect in the dirty limit as well [H, [sj • 

In the space of the parameters provided by mate- 
rial and excitation conditions the spin accumulation fea- 
tures both semiclassical and quantum properties. The 
exchange field and thus the spin accumulation can be 
engineered by electron density, excitation intensity, spin 
direction and electric currents, and should therefore be 
considered in advanced spintronic device concepts [3, fisl] . 
Our theoretical framework is general and can be extended 
to treat three-dimensional, inhomogeneous, and finite 
systems as well as the Dresselhaus SOI [40]. The elec- 
tronic structures of other non-magnetic conductors with 
significant SOI, e.g. hole gases in doped semiconductors 
or non-magnetic transition metals, are more complicated, 
but still amenable to a computational implementation of 
our method. 
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